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The dynamics of a glass transition is discussed in terms of the motion of a particle in a one dimensional correlated random 
potential. An exact calculation of the velocity V under an applied force / demonstrates a variety of dynamic regimes depending 
on the range of correlations. In a gaussian potential with correlator C(x) = a;^, we find a transition from ohmic behaviour 
(7 < 0) to creep motion V ~ exp{— const/ f^) (0 < 7 < 1). This provides a generic picture of the glass transition in systems 
where long range correlations in the effective disorder develop due to elasticity such as elastic manifolds subject to quenched 
disorder and the vortex glass transition in superconductors. 
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The driven dynamics of physical systems which can be 
modelled as elastic manifolds in quenched random media 
has received a lot of attention recently. Prominent exam- 
ples of such systems are the roughening of domain walls 
[0 , directed polymer growth |^ , motion of dislocations in 
disordered media dynamics of charge density waves 

, and surface growth in a random environment |^ . This 
recent interest was partly motivated by extensive studies 
of vortex dynamics in high-Tc superconductors. It was 
shown that this problem is related to the dynamics of an 
elastic manifold subject to quenched disorder |^-||], 

and a significant progress in qualitative understanding 
was achieved. The motion is viewed as a sequence of 
thermally activated jumps of the optimal "cell" of the 
manifold from one metastable state to the next as fa- 
vored by the applied force. In most of these systems 
the activation barriers for such a motion, which we will 
refer to as creep motion, depends on the applied force 
/ and diverges as / ^ 0, giving rise to a strongly non 
linear velocity versus applied force dependence. The un- 
limited growth of the creep barriers is taken now as the 
characteristic feature and operational definition of glassy 
dynamics. Usually the barriers grow as Usif ) f^^ 
leading to a typical velocity vs applied force dependence 
(or I-V curve) of the form V — e~^l^'^^'^^ in the creep 
regime. 

While there is now a consistent qualitative picture 
of the low temperature creep motion it is based 

mostly on scaling arguments and numerical simulations, 
and a rigorous analysis is still lacking. Although inter- 
esting new results have been obtained recently for the 
non-equilibrium dynamics of mean field models of glasses 
plf and for particle dynamics in an infinite-dimensional 
random potential |T2|Jl^ ], a general analytical derivation 
in physical models remains an unsolved problem. In the 
absence of a rigorous analysis of realistic physical 



situations one is seeking for models which are simple 
enough to be treated analytically and yet are able to 
mimic the large diversity of dynamics of real glassy sys- 
tems. A well known example is the problem of a sin- 
gle particle driven by an external force / and subject 
to a one dimensional random force field with Gaus- 
sian short range correlations. The term random force 
means that the correlator of the random potential V{x) 
Alx— w| where A is characterizes the 



is {U{x) - U{y)f 
strength of the random potential. This model is known 
as Sinai's model and has long been a subject of extensive 
studies p"5|-pl|,p[ The remarkable result obtained for this 
model is that even at finite temperature the mobility van- 
ishes below a threshold force fth ~ ^jlT . Moreover this 
system was found to exhibit anomalous diffusion, and ag- 
ing phenomena |20f| very much like what is observed in 
spin glasses p^Jl l[] . 

In this Letter we study a wider class of one dimen- 
sional random potentials with arbitrary correlations and 
find that our model can reproduce most of the existing 
regimes of glassy dynamics. The Letter is organized as 
follows. First we derive an exact formula for the veloc- 
ity in a one dimensional medium. The result is phys- 
ically transparent and simple enough to allow for ana- 
lytical treatment of very general potentials. Then we 
study Gaussian random potentials with a general corre- 
lator such that at large distances \x — y\~^ a: 



(Uix) - U{y)Y = K{x -y)^ A|x - y\' 



(1) 



-(1+7) 



and in Fourier space K{q) = U{q)U{—q) ~q^o 9 
wiih.K{z) = 2 ^^j^^{l-cosqz)k{q). We find that there 
are several regimes depending on how correlated the po- 
tential is. If correlations are short range (7 < 0), we 
recover a 'viscous fiow' (ohmic) regime where linear re- 
sponse V ^ f holds for small /. If correlations grow with 
< 7 < 1 we find a new creep regime V ^ exp(— 1//''). 
The case 7 = corresponds to the transition between 
these two regimes and we find a power law critical power 



law V versus / dependence, reminiscent of the vortex 
glass transition behaviour . Finally 7 = 1 corresponds 
to Sinai's case where V = below a threshold force. 

This suggests a deep connection with the motion of 
elastic manifold in a random medium (such as vortex 
systems in type II superconductors). In these systems 
creep behaviour arises from the interplay of the elastic- 
ity and pinning potential Q which both determine the 
creep barriers Ub- The bare pinning potential is uncor- 
related but the elasticity of the manifold generates long 
range correlations in the effective potential Ub that deter- 
mines the motion of the manifold. The role of elasticity 
as a tuning mecanism for correlations becomes transpar- 
ent upon noticing that the formation of the vortex glass is 
caused by a drastic change in elastic properties. Namely, 
the onset of shear modulus at the freezing point devel- 
ops long range correlations in the vortex system A 
phenomenological approach to describe a generic glass 
transition is to introduce a correlation length charac- 
terizing the spatial range of critical correlations, which 
diverges at the transition. The choice of K{x) on both 
sides of the transition is dictated by physical considera- 
tions. In the correlated phase (i.e glass phase) the natural 
choice is a Gaussian random potential with correlator of 
the form: 

Kix -y)^ A(( ti^)7 + log(^^) (2) 

The first term which dominates at large distance de- 
scribes long range correlations in the random potential 
and generalizes Sinai's model. The second term describes 
the behaviour at the critical point S,g — oo. Indeed 
one has 7 = and K{q) ^ l/g at the transition and 
thus K{x) ~ A In a;. The form (||) is an interpolation 
resulting from the crossover between the critical fixed 
point and the fixed point describing the glass phase. In 
the uncorrelated phase correlations are short range and 
one chooses a correlator as K{q) — 1/ \/?T~(1g)~^, i-e 
K{x) = Ko{x/£,), so as to reproduce the critical be- 
haviour for a <C a; <C Cg- An identical scenario was 
demonstrated using RG for the correlations in the free 
energy landscape at the glass transition in surface growth 
models, such as the directed polymer in d > 2 -I- 1 |14 

On a mathematical level the present model is the d = 0, 
n — I version of the problem of the dynamics of mani- 
folds of internal dimension d, in a n dimensional space. 
Remarkably, the case d — and n — s- cx) was recently 
studied by completely different techniques and seems to 
exhibit similar regimes fl^ . 

We consider the Langevin diffusion in the one dimen- 
sional quenched random potential U{x) in presence of a 
global bias / and thermal white noise r](t): 



dx{t) 
dt 



^VU{x{t)) + f + r,{t) 



(3) 



The probability density P{x,t) and the current J{x,t) 
satisfy: 



dP{x,t) 
dt 



-VJ{x,t) 



(4) 



with J[x, t) = -TS/P{x, t) + {f~ \/U{x))P{x, t). 

To derive the analytic expression for the velocity V we 
generalize to continuum models the method introduced 
by Derrida for discrete hopping problems. We con- 
sider an infinite periodic environment, i.e a periodic ran- 
dom force WU{x), of period L. The limit L — oo is taken 
at the end ||2^. One defines the periodized probability 
P{x) = J2k ^(-^ + which obeys the same equation 
as P, and corresponds to diffusion on a periodic ring 
of size L. Using the velocity for the particule on the 
infinite line can be expressed as: 



d < x{t) > 
Jt 



dxxV J 



— oo 
L 



+ 00 



dxJ(x, t) 



dxJ{x, t) 



(5) 



where J{x,t) = -TVP{x,t) + (/ -VU{x))P{x,t). At 
long time J{x, t) goes to a constant J and the asymptotic 
velocity V is exactly given hy V = JL. To find J for a 
fixed L and disorder configuration one must solve the 
stationarity equation: 



T^l^ + (V[/(x) - f)P{x) 



-J 



(6) 



with the two additional conditions P{0) — P{L) and 
dxP(x) — 1. The stationary solution with zero 
current J = 0, i.e the Gibbs distribution Po{x) = 
exp(^(— C/(x) + fx)) does not, in general, satisfy the pe- 
riodic boundary conditions. Thus V can be found from 
the solution with non-zero current: 



Fix) 



J and thus V follow from the normalization condition 
for P. In the limit L — s- oo, imposing the restriction 
C/(0) = U{L), unimportant for / > 0, (0) simplifies to: 

P{x) = - dz ^iUi.+z)~Ui.)~fz)/T 

T Jo 

and one gets the general formula for V: 
^ Jo 

valid for an arbitrary potential U{x). {A)^ denotes the 



with (7y(i)?7(t')) = 2T(5(i — i') and T is the temperature. translational average {A)^ = hm^^ooL ^ dxA{2 



The average in exists quite generally and is indepen- 
dent of the configuration of the random potential, i.e the 
velocity is self-averaging. The physical interpretation of 
in terms of Arrhenius waiting time is transparent. 
The average waiting time 1/V is a sum of Boltzman 
weights associated with the barriers the particle must 
overcome to move in the direction of the driving force. 
The highest barriers U{x + z) — U (x) with z > 0, produce 
the largest waiting times. 

The expression reveals immediately several gen- 
eral features. At large / one has V ^ f. At small 
force / 0, the response will be linear only if the 
barriers saturate, i.e do not grow at large distance. If 
the potential is uncorrelated at large distances such that 
^g(c/(x+z)-c/(x))/T^ ^ (e^/^)(e-^/^) when z ^ oo, then: 



V ^ f- 

Do-' {eU/T){e-u/T) 



(10) 



where D and Dq are the diffusion coefhcients in presence 
and in absence of disorder, respectively and the Einstein 
relation holds. When D <C Dq the V-f curve will show 
strong nonlinearity at intermediate scales where the tran- 
sition between low force and high force regime of motion 
occurs. 

We now turn to a detailed study of gaussian disorder 
with correlator K{x). Upon averaging over disorder (^ 
yields: 

K{x), 



1 1 p , , fx 
— = - c(a;exp(- — 



V T 



T 



2r2 
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(11) 



The choice of K{x) as in (2) gives rise to several regimes 
of particle dynamics depending on the range of the cor- 
relations of the random potential. 

Sinai's case 7=1: For 7 > 1 the integral in ( pi] ) di- 
verges and the velocity is zero. Sinai's model corresponds 
to 7 = 1 and appears as a marginal case where the in- 
tegral (0) diverges for f < ftf, = A/{2T^g) and V = 
while V = f — fth exactly for / > fth , in agreement with 
previous results p"5|- p^ , p9| , p0t . The system with 7 = 1 ex- 
hibits algebraic distributions of waiting times which gives 
rise to aging phenomena [20]. Therefore Sinai's model 
mimics essential aspects of the spin-glass behavior. In- 
terestingly, this V versus / dependence mimics also the 
dry friction phenomenon. 

Creep motion < 7 < 1: In the intermediate case 
< 7 < 1 one finds the " creep" dynamics regime. Defin- 
ing the dynamical exponent z ^ 2 + {A/2T'^) and the 
characteristic force fc = (A/2r2)i/7 t/^g 

one arrives 

at: 

^ ^ Hif/fc) 



y fc 



H{y) 



dv 2 exp(-?/u -I- u^) (12) 



Note that fc reduces to the threshold force fth when 7 
1. At 7 < 1 the sharp threshold disappears but at / <C 



fc the V versus / dependence shows strongly nonlinear 
behaviour with an essential singularity at small /. Using 
the steepest descent method at / ^ fc one finds: 



V = AT fc'-' (f)""^ exp(-(l - i)C-^Y) 

Ic J 



7 



1-7 



(13) 



with A ^ ^7(1 — 7)/27r. The exponential factor holds 
for any correlator behaving as a power law ~ x^ at 
large distances. The preexponential factor depends on 
details of the crossover of the correlator to the logarith- 
mic regime (x ^ ^g)- In this creep regime the linear 
response at / — + is absent and the characteristic barri- 
ers which control the dynamics diverge as 1//^. 

Critical case critical behaviour, which generically cor- 
responds to 7 = 0, can equivalently be achieved by tak- 
ing = CO in (2). The V versus / characteristics be- 
comes a power law at small force. In that case we have 
K{x) ^ A log at large x and equation (|l|) gives: 



V 



r(z - i)T2+^ 



(14) 



with z = 2 + (A/2T2). resembling the power law 
critical behaviour proposed at the vortex glass tran- 
sition. This transition separates the creep dynamics 
V ~ exp{— const/ f^) in the vortex glass state from the 
ohmic behaviour V oc f above the transition. Indeed one 
finds that the scaling function H{y) of (|lj) reproduces 
the correct scaling behaviour Hjv) ^ y'~' in the critical 
region y ^ 1 which leads to (|lj) for T/a^ 3> / '3> fc- 
Using ( p^ ) we find the critical behaviour D cx £,g^~^. 
Since in the vortex system the voltage is proportional to 
the vortex velocity and force is proportional to applied 
current the obtained behaviour of D reporduces the crit- 
ical behaviour of the resistivity p at the transition, ^g 
playing the role of the vortex glass critical length. 

Ohmic behaviour 7 < 0: For short-range correlations 
7 < 0, C/(x)2 is finite and one recovers the Einstein rela- 
tion ( p^ with D > and linear response. Let us analyze 
in detail the regime of ohmic motion, corresponding to a 
finite correlation length. We denote U{x)U{0) = A C{x), 
with C(0) = 1 and C {x) decays to zero on a length scale 



1 

V 



dzexp{-^ + ^{l-C(z))) (15) 



Two ohmic regimes, one at small force with V oc e~^^'^ f 
and one at large force with V (x f, appear. At high tem- 
perature these two regimes match smootly. At low tem- 
perature there is a sharp crossover. A depinning temper- 



ature Tp = V A separates these two different behaviours 
which we call " unpinned" and " pinned" respectively. For 
T >> Tp the particle motion is mostly unpinned and the 



smooth crossover occurs at fc ~ T/^^. For T << Tp a new 
characteristic force fp arises which marks the crossover 



between V ix e' 



/ for / << fp and F cx / for 



/ >> fp- To estimate fp we use C{z) = exp{—{z/£_y) 
which leads to fp « {T/^)y/AjT^ - \/A/^- Note that 
fp does not depend on temperature and has thus a limit 
when T 0. At T ^ the curve V versus / be- 
comes V — 9{f — fp){f — fp). It is interesting to mention 
that this behaviour can be obtained correctly using the 
perturbation theory analogous to that used for calcula- 
tion of critical currents in type II superconductors . 
Disorder induced corrections to the velocity is given by 
Sv/v = A^vT~^{j^ — arctan{j^)). Equating 6v/v fa 1 
we recover the above result for fp at low temperature 
T « Tp while for T » Tp, 5v/v never reaches 1 and 
thus fp does not exist. This analysis assume a finite 
mean squared force < (VC/)^ >~ /p^ < oo, i.e that 
the random potential is smooth. If we choose instead 
C(z) = exp(—\z / an exact calculation from (15) gives: 



V = 



T 



C7(/$/T,l)'2r2 



where 7 is the incomplete Gamma function. In that case 
the average force is infinite and at T = the particle 
remains pinned at any applied force. 

The above model may apply directly to a string (di- 
rected polymer) of length L in d dimensions driven by its 
tip at z = over point like impurities. An example is a 
flux line in a superconductor in presence of an external 
current J. Because of screening at small fields the exter- 
nal Lorentz force of modulus / oc J is exerted only on the 
end of the line, over a length of order A. Let F{r) be the 
free energy corresponding to equilibrium positions r of 
the tip at z = with the other end z = —L fixed at r = 
for large L. One can view the tip as a particle driven in 
the rough potential F{r) with {F{r) — F{0)y ~ and 
7 — 20 /(^. Moving the tip by r provokes a reorga- 
nization of the string up to length z = r^/'» away from 
the tip. It thus leads to correlated barriers and to the 
V-f characteristics found above. Since C ~ 6/(8 + d) 
and 6* = 2C - 1 one has 7 « (4 - d)/3. Thus 7 = 1 in 
d = 1 -|- 1, an exact result, and 7 « 2/3 in d = 2 + 1. If 
this picture is qualitatively correct a flux line driven by 
its tip in d = 1 -|- 1 would experience a threshold. Below 
fc the motion of the tip would be as r ~ t^^^" and V = 
and V ^ exp(— 1//2) in d = 2 + 1. Of course the true 
motion is more complicated and the string might not be 
assumed to be in equilibrium. 

Note that in one dimensional motion the velocity can 
be dominated by rare large barriers, as is the case for 
the exponent fi in (13), whereas diffusion at zero force 
/ = 0, given ||,|5fby x - (logt)^/^, is dictated by 



typical barrier Eh 



In higher dimensional space. 



such as the configurational space of a string, there are 
many paths in parallel from one point to another which 



may allow to avoid atypically high barriers. This could 
effectively cut the tails of the distribution of barriers, 
and the gaussian assumption may not be justified. This 
effect can be accounted for within our one dimensional 
model by allowing more general distributions for barriers. 
Averaging (|^) over disorder, one can write: 



1 

V 



dz 



J dEbPiEb,z)exp{{Ei 



\-fz)/T) (16) 



where P(Ei,,z) is the probability of a barrier Ef, = 
U{x + z) — U{x) corresponding to separation z. The pre- 
vious calculation corresponds to a Gaussian P{Eh,z) ^ 
exp{—Eb^/2K{z)). A more general form is: 



P{Et,z)^z-'>^^Q{Eb/z'>^^) 



(17) 



The behaviour at small / is governed by the large barri- 
ers. Thus we suppose that Q{u) ^ exp{—u^/2A) at large 
u. Then using the saddle point method one obtains: 



y ~exp(-cr-(i+'')/"^) 



7(5 



1 (18) 



withiy= l/((5-l-^) andC = i/-i(2/7'5)^+''(2A)-i/>^. 
This is the most general case containing the Gaussian 
case for S — 2. For 6 00, i.e when the distribution has 
no tails, the typical barriers dominate and n = iityp = 
7/(2 — 7). This calculation shows however that even for 
rapidly falling tails, i5 large but finite, the exponent fx can 
be different from the value ^typ obtained by considering 
typical barriers. 

In conclusion, we found stretched exponential veloc- 
ity versus force characteristics V ~ e:xjp{— const/ f^) in 
a simple model of a particle diffusing in a ID random 
potential with correlated barriers. There is a transition 
between this creep motion and ohmic motion, character- 
ized by power law critical behaviour at the transition. 
This transition occurs as the correlations in the random 
potential change from short range to long range. We 
stressed the physical analogies with the creep motion of 
driven flux lines in a glassy state, also characterized by 
correlated barriers. The transition found here is similar 
to the vortex glass transition which can thus be viewed 
as a transition from the short range (in the liquid) to 
the long range correlated behaviour (in a glassy state) 
of the effective potential seen by a vortex configuration. 
The above results also establish a connection between 
spin and vortex glasses. Our model exhibits both creep 
behavior specific to vortex glasses and aging, which is 
recognized as the essential characteristic of spin glass dy- 
namics, upon tuning the degree of correlation of random 
field. 
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Note added: After submission we received a preprint by 
S. Scheidl who also derived Eq. 8. in a different context. 
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